Alternative approach to computing transport coefficients: application to conductivity 
and Hall coefficient of hydrogenated amorphous silicon 
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We introduce a theoretical framework for computing transport coefficients for complex materi- 
als. As a first example, we resolve long-standing inconsistencies between experiment and theory 
pertaining to the conductivity and Hall mobility for amorphous silicon and show that the Hall 
sign anomaly is a consequence of localized states. Next, we compute the AC conductivity of amor- 
phous polyanaline. The formalism is applicable to complex materials involving defects and band-tail 
states originating from static topological disorder and extended states. The method may be readily 
integrated with current ab initio methods. 



The atomistic understanding of the electrical 
conductivity [H and Hall coefficient 0-0] are key 
unsolved problem in the physics of amorphous 
semiconductors iS]. The challenge is due to two 
factors: (1) a complex array of localized states of varying 
physical origin; and (2) at moderate temperature, 
both localized states and extended states are accessible 
excited states [9l.[lo|. 



Previous work on the conductivity [UJ, [12] and Hall 
coefficient of amorphous semiconductors has in- 

volved either Miller- Abrahams theory or small polaron 
models'^. Neither approach is ideal for amorphous semi- 
conductors, with their complex menagerie of localized 
states due to defects, and tail states due to topological [18 
or chemical disorder Q and electron-phonon couplings 
strongly dependent upon carrier localization 

In this Letter, we first develop a novel formalism for 
the linear response to a mechanical perturbation's^. The 
method properly includes the four possible transitions be- 
tween extended or localized initial and final states. The 
full results and many intermediate details are provided 
in Ref. 2J|; here we require only transitions between lo- 
calized states. Our work resolves the puzzle of the sign 
anomaly of the Hall mobility in a-Si:H:, provides the tem- 
perature dependence of the transport coefficients, and ac- 
curately predicts the ac conductivity of polyaniline. The 
new formalism is expected to be applicable well beyond 
the examples of this Letter and lends itself to inclusion 
in current ab initio schemes. 

In the small polaron hopping regime, the Kubo linear 
response formula [20j has been used to compute conduc- 
tivity and Hall mobility l^, Il5l4l7| . The key mathemat- 
ical obstacle to computing the current-current correla- 
tion function is the imaginary time integral, 12, 15-l3| 
which results from the commutator between the micro- 
scopic current and density matrix poj. For a "mechanical 
perturbation" (for which an external disturbance may be 
expressed with additional terms in the Hamiltonian[20|), 
we can avoid this troublesome commutator, and within 
this picture the local density operator p of a quantity 
(charge, energy etc.) can be easily constructed[22|. 



the mechanical perturbation, from which the microscopic 
local density p{r,t) = {^'{t)\p\'^'{t)) is obtained. Next 
we calculate dp{r,t)/dt by means of the time-dependent 
Schrodinger equation ihd'^'{t)/dt = H'{t)^'(t), where 
H'{t) is the total Hamiltonian of [system + mechanical 
perturbation]. The third step is to apply the local den- 
sity (charge density, energy density) continuity equation: 
dp{r,t)/dt + \7- jm(r, t) = 0; the microscopic response 
jm (current density, energy flux etc.) is then expressed in 
terms of ^''(i). Practically speaking, ^''(i) may be com- 
puted to the required order with perturbation theory. By 
substituting ^'(i) into the expression for j„i(r, t), one can 
obtain the microscopic response to the required order of 
mechanical disturbance. Spatial and ensemble average 
are taken at the final stage. The desired transport coef- 
ficients can then be extracted from the ensemble average 
of the spatially averaged fiux j. 

Since the state of the system is completely determined 
from the initial conditions, averaging over initial state 
can be delayed until the final stage. Thus we can avoid 
the commutator between flux and density matrix, i.e. 
the integral over imaginary time. 

Consider then, a system with electrons and TV 
nuclei in the presence of an electromagnetic field with 
potentials (A, 0), the charge density of state 4*' at r 
p'{r,t) — J dT'i>'*p{r)'i>' , where the arguments of wave 
function Vf' are (ri • • • r^r^ ; Wi • • • Wjv'; t), Wi is the 
position of the first nucleus, etc. p{r) = J^j ^^(.^ ~ 
Tj) — J2l ZLeS{r — Wl) is the charge density operator. 
dr = dri ■ ■ ■ dr^^dWi ■ ■ ■ dW j\f is the volume element in 
configuration space. The evolution of the state is deter- 
mined by the time-dependent Schrodinger equation for 
which the total Hamiltonian includes the interaction be- 
tween system and external field. The contribution from 
the electrons is: 



A(r,i) J dr'*'**'. 



We first average p over a state ^''(t) of the system with where the arguments 



of 



(1) 



are 



2 



(r,r2,--- ,r7vjWi,W2,--- ,WA/-;t), dr' = 
dr2 • ■ • drAr^dWi(iW2 • • • dWN„ ■ In Eq.((l|), the anti- 
symmetry of 5*' under exchange of particles was 
used, and the Coulomb gauge V • A(r, t) = was 
adopted to simplify the expression. j„i(r, i) is gauge 



invariant [l4 |23|. Without the nuclear coordinates, 
Eq.(IT|) is identical to the form used by BCS to compute 
paramagnetic and diamagnetic currents for supercon- 
ductors in which Bloch states are not radically modified 
by the electron-phonon (e-ph) interaction [i^. Eq.([T]) is 
a generalization to arbitrary strength of e-ph interaction, 
and may be used for the localized carriers in amorphous 
semiconductors or the polarons in ionic and molecular 
crystals. 

We now apply Eq. ([T]) to compute the conductivity and 
Hall mobility of an amorphous semiconductor. The ratio 
of the second term to the first term is eA/p (p is elec- 
tron momentum), so that the contribution from the 2nd 
term may be neglected. Since the carrier concentration 
is low in lightly-doped amorphous semiconductors, one 
can invoke the single-electron approximation, and ^'(t) 
for jm may be replaced by the single-electron wave func- 
tion 'ip'{r;xi,X2,--- , X3j\f; t), where r is the coordinate of 
the carrier, xi,X2, - ■ ■ ,X3j\f are the displacements of 3Af 
vibrational degrees of freedom. 

Using perturbation theory, one can expand ip'it) to the 
required order of external field 13, 21 1. We use A with 
subscripts to label localized states, denote the coupling 
between two localized states (f)A and (fiAi caused by ex- 
ternal field as Ja^a ~ J dr(t)*^_^hfm4'A, where hf„i — 
(ihe/m) A{r) ■ Vr + e'^ A'^ {r) / {2m) + e(j){r). The spatially 
averaged microscopic current density to second order of 

jfield -g 



j(s,i) = - 



Nehe 
mf2s 



dr 



3Af 

([[dXj 



+ Im(V'(°' V?/'^^'* - V^^^* VV-'"^ + V^^' VV'^^'*)}, (2) 

where t/j^-^^ is change in state to order jl^^'-'^^ and ■0^^^ 
is change in state to order [J^""^'"*]^, where ip'^'^^t) is the 
state of carrier at time t without external field, ilg is the 
"physical infinitesimal" volume of Kubo^24|. Because the 
initial state of the phonon-dressed carrier is unknown, 
we need to average jkik = x,y,z) over initial phonon 
distribution and single-electron states. 

To compute the conductivity, we only require the order 
J field f^Qj-frj^ of Eq.Q. If one applies a voltage drop across 
the material, the potentials are A = and (f> = — 2Eo • 
rcoswt. Using perturbation theory, one can compute 
■0^°^ and V'*-^'' to order J^. Substitute ■i/'^^^ and ip'^^^ into 
the first term of Eq.([5]), the conductivity from the LL 



transition is 1211: 



Im f7Q/3(a;) 



lmi^[lAiA+ ± IaiA- 



AAi 



x/(£;0)[l - f{EX)]vXA{El - E"aX'«a. <a), 



^^Im ^ 



AA1A3 



f{E'Am~f{EX)][l-f{E'iJ] 



'^{Ea ^ Ej^J h Ja3a[Ia3AiA+ ± IA3A1J 



X -^ma3>a\a 



(3) 



A2A1A 



{Q1A2A1A+ ± QlA2AiA-)f{EA)[^ - I{E% 



'2hn^ 



J2 lH<A2-vl2A)Jl2MiEl-E'AX'iyMAr 
A2A1A 



{Q2A2A,A+ ± Q2A2A,A^)f{E'A)[l - /(i^^)] 

Re(jjfe(w) (j, fc — x,y,z) takes the positive sign and 
Imcrj7j(w) takes the negative sign. Unlike previous 
theories [9[, the dc conductivity may be directly extracted 
from Eq.Q without a limiting process. Here, v\_^j^ = 
/ dv(t>\^{pk/m)4)Aa.^dw\j^^ = / dY(l)A{pj /'m)(l)\^ arethe 
velocity matrix elements of a carrier. f{Ej^) is the Fermi 

distribution. Ja^a = I dr(l)\^J2p^DA^^'^ ^ ^p")^^ 
transfer integral from c/ja to (j)Ai- In Eq.([3]), the first 
term is order J*^ contribution, the second term is order 
contribution. Here, IaiA±{'^,T) has dimension of time 
(denote its order as t^^^), and reflects the time evolution 
of ■(/''(i) in a field-driven 2-site transition, cf. Eq.(Cl) in 
[2H . The characteristic time t^ ^ may be interpreted as 
the mean free carrier time. Similarly the two-fold time 
integrals Ia3AiA±, QiA2AtA±, and Qia2AiA±, have di- 
mension of [time]^, reflecting the time evolution of ^p'{t) 
in a 3-site transition induced by transfer integral J; they 
have same order of mag nitude [i^ cf. Eqs.(C3,C7- 

C9) in 21]. JA3A[tA3AiA\'^ the mean free time in 
a 3-site process. The interference between jf^^^'^ and J 
is displayed between two components (reached through 
different paths) of the final statelHj. At high tempera- 
ture kgT > tiOj (uJ is first peak in phonon spectrum), the 
time integrals Is and Qs can be approximately computed 
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FIG. 1: AC conductivity of polyaniline as a function of 
frequency at T=300K: star symbols denote experimental 
data[2^, solid line is from the first three terms in Eq.(|4]). 



using the method of steepest descent (see Appendix D in 
211). For example: 

/9'i("AiAT'^)/2-/3AAiA/4 



iAiA± — e 



{C 



n=0 



2 



cjc«cosh— — ] }, 



(4) 



where 9^^ is the shift in origin of the a*'' normal mode 
caused by a carrier in state (t>A, , Xa, a = \ X^a f^a{Oa^ — 
e^f and C = hr'^kBTXA^Am, lill- The mean free 
time decreases with increasing XaiA and energy differ- 



ence: 
The 



average 



K^bTXaiaY 
mobility /i 



l/2g-/3?i(i^:ij^T'^)/2-^AAiA/4^ 

is defined by a^j. = 



e^(A^e/f^s)M- From Eq.®, one can see that /i de- 
pends on the energy distribution and spatial distri- 
bution of localized states. A typical value of ji 



can be estimated: 



"AiA 



tAiA {Ea 



A 
-1 



^Ai 



VA,AVA^AAtAsA,Afl ^Ja,a) {E% - E%J . 

As a test, we apply Eqs. pi4p to the frequency 
dependence of the ac conductivity in polyaniline at 
T=300K 25]. The Austin-Mott uj°-^ lawfl does not ac- 
curately fit experiments [2^. In Figll] we fit the datapsj 
with the first three terms (a quadratic polynomial of uj) 
in Eq.dH). Because a;,„ax = 10'^IIz<< T = 300K, factors 
^±0huj ~ ]^ (Jo not play a role in the low frequency regime 
hio « ksT. 

To describe the Hall effect, one needs ip'{t) to second 
order of one j/*'^''' includes electric field E^, and 

another includes magnetic field B^- After substituting 
^(0)^ ^(1) and V'^^^ into the second term in Eq.Q, and 
averaging over various initial conditions, <7yx is deter- 



mined from jy, see Eq.(49) in [21|. The primary tem- 



perature dependence of ayz is included in the time in- 
tegrals which are obtained from integrating out vibra- 
tional states. The 3-site processes result to two-fold time 



integrals 2l| which have dimension [time]^. The time in- 



tegrals have the same order of magnitude Saa2Ai^ where 
SAA2A1 is the characteristic time for the 3-site processes. 
Similarly the 4-site processes result to three-fold time 
integrals [2l|, the order of magnitude is s\a3A2Ai' where 
SAA3A2A1 is the characteristic time for the 3-site pro- 
cesses. SAA2A1 (3-site processes) and [fi-~^ Jsaa3A2AiV^'^ 
(4-site processes) may be explained as the mean free 
times in presence of magnetic field. Applying fields Ey 
and Bz, (Jxy is inferred from jx- The order contri- 
butions come from various 4-site processes, their order 

is Uee^ {h-^JA3ASAA3A2Ai)f^~'^ [BzeL^^^^ /m]XAiAVA3A2 , 

where Lf^^^ = J dr(j)\ Lz<j)Ai is the matrix element of 
the z component of electronic orbital angular momentum. 
The widely used 3-site process in the literature 0, 0, [13- 
[l6| is a special case when Ai = Aj,. The order J° 
contributions comes from 3-site processes. Their or- 



der is 



h [B zcL z I Tn)xv . It roughly cor- 



responds to the 'interference' contribution (involving 
two sites') IIGII . They are neglected in some other 
calculations |l4 [l7j . The ratio of the order terms to 

jO terms is {^~^ J A3As\a3A2A^s~aA2A^ ~ 1- 

Amorphous semiconductors are isotropic, so that one 
may estimate the Hall mobility as: [in = B^^Gxy/oxx ~ 



j^XVSaA2A 

For a-Si:H 



AiA^ where 0:^3^1 = 
^ and RaiA ~ 5 



ihv / 



i{E\-Ei^y 



A3A\-^A 

lOA, J ~ 0.02eV, one 
The temperature 
dependence of [lr may be obtained from those of Uxy 
and axx'- 



has [iH ~ 0.1 — 0.2cm^Volt ^scc 



jXVSaA2Ai^AiA 



exp{- 



5A1 



r,A\2 



Y)- (5) 



At low frequency, the phonon spectral density oc o;^, so 
that the sum in Eq.(l5]) converges. Figl2] gives the Hall 
mobility vs. temperature for n-type a-Si:H. Both Eq.([5]) 
and the Friedman- Holstein result [lij. are roughly con- 
sistent with experimental data[8|. E^ a — Xaia{^ + 
Aia/ XaiaY / is estimated from typical parameters 
AGaiA = O.OSeV and Xa,a = 0.2eV for a-Si^. 

We may also demonstrate the undetermined sign of 
the Hall voltage from Ehrenfest's theorem. The expected 
value for the acceleration of a carrier is: 



-/.rr(t)(-)^'(t).-|-. 



J dTtlj'*i}j'V X B 



+ dTi,'*{t)[q^ - V,V{v, {W„})]V'(t) (6) 

~ [ dTiP'*B X (-iWr?/-') + / dTi;'*-^(B X AU', 
J J 



4 



■0.09 




potential to be implemented with current single-particle 
ab initio simulations, and requires only the eigenvalues 
and eigenvectors of the single-electron Hamiltonian, dy- 
namical matrix, and quantities easily derived from these. 

We thank the Army Research Ofhce for support under 
MURI W91NF-06-2-0026. 



FIG. 2: Hall mobility vs. temperature: square symbols for n- 
type a-Si: dashed line from best fit of Friedman- Holstein 
formula 14], solid line is from Eq.((5]). 



where q is the charge of the carrier. If a system has only 
extended states, because the mean free time of a carrier 
is much larger than the effective interaction time with 
phonons and defects, the drift velocity is along the di- 
rection of gE for a nearly free carrier. The direction of 
the average magnetic force in an extended state (diag- 
onal element) is the same as the classical one gE x B, 
and the sign of the Hall voltage is normal. Because the 
force exerted on a carrier by the external E field is much 
weaker than the binding force of the disorder potential, 
the time average of 'm~^ J <f>\_^ {—ih\7r — qA)(j)Ai in a lo- 
calized states (pAi is zero: no magnetic force acts on a 
localized carrier. The mean trapping time that a carrier 
spends in a localized state before making a transition to 
other states is[i3 ~ hJ-^(\kBTY'^e^-'^^'^ (high tem- 
perature) or ~ hJ~^{AG) (low temperature), where AG 
is the typical energy difference between the final and the 
initial states, and A is the reorganization energy for a 
transition. The mean transition time needed for a transi- 
tion event is ~ mdR/h, where d is a typical bond length, 
R is the distance between two localized states for LL 
transition, R is localization length for localized-extended 
transition. The mean transition time is much shorter 
than the mean trapping time in a typical localized state. 
Comparing with the transition speed h/{md), the speed 
qEdR/h obtained from external electric field E during 
the transition time is negligible. The magnetic force suf- 
fered by a carrier during a transition is along the direction 
of X B, where v is the transition velocity of the car- 
rier. V does not have any relation to the direction of E. 
If one applies E field along the x axis and B field along 
the z axis, a direction dependent Hall voltage should be 
detectable along any direction in the yz plane, not only 
along the y axis. Checking this prediction would be a 
critical test for this work. A recent experiment [27] shows 
that the signs of Hall voltage in several a-Si:H films are 
not always reverse to those expected from gE x B, the 
present analysis seems to agree with this observation. 

We show that localized states are responsible for the 
anomalous sign of the Hall voltage. The method has the 
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